Ion-water clusters, bulk medium effects, and ion hydration 
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Thermochemistry of gas-phase ion-water clusters together with estimates of the hydration free 
energy of the clusters and the water ligands are used to calculate the hydration free energy of the ion. 
Often the hydration calculations use a continuum model of the solvent. The primitive quasichemical 
approximation to the quasichemical theory provides a transparent framework to anchor such efforts. 
Here we evaluate the approximations inherent in the primitive quasichemical approach and elucidate 
the different roles of the bulk medium. We find that the bulk medium can stabilize configurations 
of the cluster that are usually not observed in the gas phase, while also simultaneously lowering 
the excess chemical potential of the ion. This effect is more pronounced for soft ions. Since the 
coordination number that minimizes the excess chemical potential of the ion is identified as the 
optimal or most probable coordination number, for such soft ions, the optimum cluster size and the 
hydration thermodynamics obtained without account of the bulk medium on the ion- water clustering 
reaction can be different from those observed in simulations of the aqueous ion. The ideas presented 
in this work are expected to be relevant to experimental studies that translate thermochemistry of 
ion- water clusters to the thermodynamics of the hydrated ion and to evolving theoretical approaches 
that combine high-level calculations on clusters with coarse-grained models of the medium. 
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I. INTRODUCTION 

Ion-water clusters represent the transition between an 
ion in the gas phase and a fully hydrated ion. Hence 
these systems have been extensively studied to under- 
stand the thermochemistry of hydration and to infer 
how bulk properties emerge from properties of ions in 
clusterai"— . Most notably, attempts have been made to 
infer the hydration free energy of an isolated ion, a ther- 
modynamic descriptor of the nonideal interactions be- 
tween the ion and water, by analyzing the thermochem- 
istry of forming ion- water clusters X[H20]if2r2i, where X 
is the ion. 

Thermochemical data alone is inadequate for inferring 
the structural characteristics of hydration. In this regard, 
theoretical calculations of cluster formation have proven 
useful. (For example, see Refs. [l^ - [l9l .) Moreover, such 
calculations together with hydration free energy calcula- 
tions of the clusters and water ligands have been used 
for estimating the hydration free energy of single ions. 
In calculating the hydration free energies, almost always 
a continuum model for the bulk solvent is assumed and 
either onei^ or an ensemble^'' of gas phase configurations 
of the cluster is used for obtaining the hydration free en- 
ergy of the cluster. In some studies the gas-phase cluster 
is also allowed to relax in response to the model bulk 
medium^^ . 

An important development in the statistical me- 
chanics of hydration has been the quasichemical 
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of the potential distribution theorem. 



This approach provides a rigorous, statistical mechani- 
cal framework to relate the thermodynamics of solute- 
water clusters to the bulk hydration of the solute. In 
practical implementations, approximations to the quasi- 
chemical theory are invariably made. In the primitive 



quasichemical approximation, — the adjective primitive 
indicates neglecting the role of the bulk medium on the 
local clustering reaction — the thermochemistry of ion- 
water clustering in the ideal gas phase, obtained using 
standard quantum chemistry approaches, is coupled with 
estimates of the response of the bulk medium, treated ei- 
ther as a dielectric continuuni^^"— or as a discrete molec- 
ular solvent;^. 

Within primitive quasichemical theory, the hydration 
free energy is estimated by varying the number of water 
molecules in the cluster and finding the optimal coordina- 
tion number that minimizes the hydration free energy of 
the ion. This cluster variation approach has proven suc- 
cessful in establishing the hydration structure of small, 
hard ions such as Li+(aq)2^, Na"*" (aq)2Si, and Be^"'"(aq)^; 
for these cases, the optimal coordination structure pre- 
dicted by primitive quasichemical theory is in good agree- 
ment with the most probable coordination observed in ab 
initio molecular dynamics simulations (AIMD). Similar 
agreement is also seen for Mg^"*" , Ca^"^ , and some of the 
transition metal dicationa^. Given the then limitations 
of small systems, short simulation times, and uncertain 
quality of the underlying electro density functionals in 
AIMD simulations (for example, Refs.i33andi34i), primi- 
tive quasichemical theory proved useful in cross-checking 
the simulation results itself. (The predictions for H+ — 
and HO^"^^ are also in fair agreement with AIMD simula- 
tions, but consensus remains elusive^li^. These systems 
also challenge both theory and simulations because of the 
need to describe nuclear quantum effects.) 

For the soft K+ ion, primitive quasichemical pre- 
dicts an optimum coordination with four (4) water 
molecules^. But the results from AIMD simula- 
tions are less conclusive. Some studies identify an 
inner-coordination number of four and a second outer- 
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population with two additional water molecules within 
the nominal first hydration shell of the ion'^'^''^^, but oth- 
ers make no such distinction^. Beside K+, it has also 
been well appreciated that for some of the halides, op- 
timum gas-phase clusters can show intcrmolecular bond- 
ing between the coordinating shell water molecules^^^, 
a feature that is not usually observed in the coordina- 
tion structure for the hydrated ion in the liquid. Thus 
the primitive quasichemical approach for these ions has 
not been successful in reproducing the hydration struc- 
ture and thermodynamics. The limitations in describing 
the hydration of soft cations and anions suggests that 
the bulk medium plays an important role in the hydra- 
tion structure of these ions. Developing the framework 
to understand this effect is the objective of this article. 

In Sectionini we present the quasichemical theory, elu- 
cidate the role of the medium, and highlight the physical 
consequences of the primitive quasichemical approxima- 
tion. In our previous work**^ we showed the importance 
of occupancy number variations of water in an empty 
coordination sphere for understanding an ion's coordina- 
tion structure. Developing those ideas further, here we 
find that the bulk medium promotes a better packing of 
solvent molecules around the ion, leading to a decrease 
in the contribution to the hydration free energy due to 
interactions between the ion and the solvent molecules 
within the coordination sphere. The medium stabilizes 
configurations which are otherwise not observed in iso- 
lated ion-solvent clusters. Without proper account of 
this effect, the predicted optimal coordination state is 
typically lower than the optimal coordination observed 
in simulations, and appreciating this effect is necessary 
in calculating hydration thermodynamics as well. 

II. THEORY 

A. Quasichemical theory 

We present only the main elements of the quasi- 
chemical theory; more extensive discussions are avail- 
able elsewhere2i"— . We define the coordination sphere 
of radius A around the ion, a, and by so doing, sepa- 
rate the local, chemically intricate ion interactions with 
water molecules (the solvent ligands) within the coordi- 
nation sphere from the interaction of the ion with the 
medium (the bulk) outside the coordination shell. The 
chemical potential can then be written as^i^— 

/^/i^ =lnxo-lnpo + /3//^"(0), (1) 

where j3 = l/fcsT, fc^ is the Boltzmann constant and 
T is the temperature, xq is the probability of finding 
zero (0) solvent ligands within A. (The position of the 
water molecule is defined by the position of the oxygen 
atom. The ion in the n*"^ hydration state has n-ligands 
within A.) The contribution to the free energy due to 
ion-bulk interactions when the ion is in the n = state 
is fJ.'^in = 0) = /i°'^(0). The probability of finding zero 



(0) solvent ligands within A when the ion-water interac- 
tions are turned off is pq. This factor accounts for the 
free energy of creating an empty coordination sphere in 
the solvent and thus accounts for the packing (steric) 
contributions to fj,"^. 

The quasichemical form is obtained by noting that xq 
is specified by chemical equilibria of the form^ir— i^i^ 

a[H20]o(aq) + nH20(aq) ^ a[H20]„(aq) , (2) 

where a[H2 0]„ denotes the ion plus n- water cluster 
within the coordination sphere and 'aq' indicates that the 
clustering reaction is in the presence of the bulk medium. 
Given the equilibrium constant of the above reaction is 
Kn and the density of water molecules is p^, we imme- 
diately obtaiit2i-24d2i43 

\nxo = -H^ + J2K„pl). (3) 

n>l 

Note that the excess chemical potential of the cluster 
a[H2 0]o is just 

/3A^a"»o =/?Mo'(0)-Wo. (4) 

(For ease in writing, we will denote II2O by w in all 
subscripts.) 

Approximating Eq. [3] by its maximum term, we get 

«-lni^nP:+/3/i^"^„. (5) 

Kn is given by the ratio of partition functions of 
the cluster to individual molecules and can be challeng- 
ing to calculate with full account of the surrounding 
mediu m^^i^^i^^i*^^" — . But the equilibrium constant Kn ^ 
of the clustering reaction 

a[H20]o(g) + nH20(g) ^ a[H20]„(g) , (6) 

in the ideal gas phase is more amenable to theoretical cal- 
culations. In Eq. [SI 'g' indicates that the solvent outside 
the coordination shell is non-interacting (or ideal); that 
is, the reaction is performed in the ideal gas phase. Since 
the free energy change for Eqs.[2]and|6l with appropriate 
choice of standard concentrations, are just —k^T In Kn 
and —kBTliiKn \ respectively, we find that 

ksT In -% = p'^^.^„ - pZ^^ + n/.- . (7) 

Kn 

The first two terms on the right hand side represent the 
hydration free energies of the clusters containing waters 
and n waters, respectively, and ia!^ is the hydration free 
energy of a water molecule. 

A similar development can be pursued when the coor- 
dination shell is empty, which is equivalent to the case 
when the ion and solvent do not interact. Denoting the 
noninteracting solute by •, the analog of Eq. [7]is 

fcsT In 4^ = pl^^^ - pl^^^ + npZ^ . (8) 
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Note that po is related to Kn~^ exactly like Xq is related 
to Kn (Eq. [3]), and equilibria analogous to Eqs. [2] andlB] 

specify Kn and Kn \ respectively. 
From Eqs.[5]and [7] we obtain: 



Mo 



-kr,T In K^°h 



(9) 



Eq [9] not only permits the calculation of but it also 
suggests an approach to identify the optimal coordination 
state. For different n, we can compute Kn^; often highly 
descriptive forcefields, including ab initio potentials, are 
used to model the thermochemistry of the gas-phase reac- 
tion. The presence of the bulk medium outside the coor- 
dination shell is then corrected a posteriori by adding the 
free energy of hydrating the cluster and subtracting the 
free energy of transferring the requisite number of water 
molecules to the gas phase from the liquid. A continuum 
or dielectric model of the solvent is often assumed for 
this purpose. The coordination state that minimizes 
is then identified as the optimal coordination state^^"— . 
In the primitive quasichemical approximation, the con- 
figuration of the gas phase cluster is not allowed to relax 
in computing the hydration free energy. To better un- 
derstand this approximation, we first identify the role of 
the bulk medium on local ion-water clustering. 



B. Role of the bulk medium on local ion- water 
clustering 

The ratio of equilibrium constants Kn and A'„ was 
earlier— shown to be 



K„ 



(10) 



where Wn is the free energy of forming the ion-water 
cluster in the presence of the bulk and /i^^ (n) is hydra- 
tion free energy of the ion in its n*^ coordination state. 
Specifically, 



-/3Mr(n) 



-fie 



(11) 



where e is the interaction energy of the ion with the rest 
of the medium, (. . .)o indicates averaging when the ion 
and the medium are thermally uncoupled (denoted by 
the subscript 0), and (. . . |n)o indicates that only cases 
with exactly n solvent molecules within the coordination 
shell are considered in averaging. 

By parsing the interaction energy e into a local piece, 
^locai, obtained by considering ion-interactions with the 
n coordinating water molecules, and the remaining long- 
range piece, eir, we have 



(e-/3^iocai|n)o 



-/3eio 



|n)o • (e-'5""» local 



(12) 



where we have made use of the rule of averages^^— , 
in rewriting the second factor on the right; (. . . |n)iocai 



indicates averaging such that the n-coordination shell 
molecules are thermally coupled with the ion, and the 
bulk medium is uncoupled from the ion but is coupled 
with the ri- water molecules. 

Identifying the first factor on the right in Eq. [T^] by 
g-^?aq(n) second factor by e~'^^outer(")^ we have 



,(")p-/3M=;to,(«) 



(13) 



We emphasize that e"^^'"'*-"-' accounts for the contribu- 
tion to the chemical potential due to the interaction of 
the ion with solvent ligands within the coordination shell 
in the presence of the medium outside^ and e"''''""*""'^"^ 
accounts for the contribution due to the interaction of 
the bulk medium with the ion when there are n solvent 
ligands present inside the coordination shell. Note that 
when n = 0, /ioutcr(O) = Mo"(0) (Eq.©. 
Thus we find that 



K. 



Kn 



(14) 



(Our earlier work^ had an error in factoring Kn/Kn- 
The error is discussed in Appendix 1X1 for completeness.) 

Now if the medium outside the coordination shell was 
an ideal gas, 



K 



(0) 



K 



^ g-/3«g(") 



(15) 



From Eq. [14] and [151 '^o get 



Kn/Kj- 
Kn/K'i 



(0) 



(0) 



= p-/3(A'SStcr(«)-/'r(0))p-/3(«a<,(n)-«g(n)) 



(16) 



and this together with Eqs.[71|Sl and [TBI gives 

M^"^„ = Uin) ~ + + Mr.„ (17) 

Finally, substituting in Eq. [SI and noting Eq. SI we get 
= -kBTlnKl^^pl + Caq(n) - Cg(n) 



C. Primitive quasichemical aproximation 



(18) 



EquationlTTlclearlv identifies three different roles of the 
solvent medium and also helps us induce the physical con- 
sequences of the primitive quasichemical approach. The 
quantity ^aq('T-) — Cg("-) is the change in the free energy 
due to the local ion solvent interaction upon coupling the 
gas-phase cluster with the bulk medium. It accounts for 
the energetics associated with the relaxation of the clus- 
ter in the presence of the medium. Utilizing gas phase 
geometries completely ignores this effect. Since the con- 
figuration of clusters with a high coordination state can 
be expected to be more sensitive to the presence of the 
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bulk medium than configurations with lower coordina- 
tion states, ^aq(f^) — '?g("-) is expected to be larger for 
higher coordination states than for lower ones. 

The quantity /^^^^^^(n) accounts for the interaction of 
the bulk medium with the ion when there are n solvent 
ligands inside the coordination shell. Since the config- 
uration of the cluster will change in the presence of the 
bulk medium, this quantity is also inadequately described 
when we do not allow the gas-phase cluster to relax. We 
do expect, however, that since only long range ion-solvent 
interactions contribute to //^^^^^(n), this quantity can be 
well captured by dielectric models of hydration. 

The final term is the hydration free energy of an n- 
water cluster without the ion. It consists of packing inter- 
actions which are also typically ignored in the primitive 
quasichemical approach. 



III. METHODS 

Liquid water and ion water systems are studied 
under NVT conditions using Metropolis Monte Carlo 
simulations^!^. The cubic simulation cell comprises 306 
water molecules for the pure water system; ion- water sys- 
tem consists of an additional ion which is held fixed at 
the center. The box volume (L^) is adjusted such that 
the number density in each case is 33.33 nm~^. Wa- 
ter is modeled with the SPC/E potential''^. Ion pa- 
rameters and the magnitude of coordination radii are 
taken from our earlier study Electrostatic interac- 
tions are modeled by the generalized reaction field (GRF) 
approacb^"-^. Both Lennard- Jones and electrostatic in- 
teractions were truncated at L/2. 

As in our earlier work^^, Na"*", K+, and F~ ions were 
studied. Since the trends for F~ follow those of the 
smaller Na+ ion, for clarity we only present results com- 
paring Na+ and K"*". The parameters for the ions were 
obtained from Ref. (soj. 

Simulations are carried out for 6x10^ sweeps. Each 
sweep consists of one attempted translation or rotation 
of each water molecule. The first 3x10^ are set aside for 
equilibration. During the equilibration phase the max- 
imum angular deflection and linear displacement of the 
water molecule is optimized to yield an acceptance ra- 
tio of 0.3. The optimized values are then held fixed for 
the next 3x10^ sweeps. Configurations of the system are 
saved every 10 sweeps for analysis. 

A similar simulation strategy was used for simulating 
ion- water and neat water clusters in the absence of the 
bulk medium. The solvent ligands comprising the clus- 
ter are restricted to a sphere of radius A. Clusters were 
simulated for 6x10^ sweeps with the first 3x10'' sweeps 
used for equilibration and the next 3x 10^ sweeps used for 
production. During the production phase configurations 
of the clusters were saved every 10 sweeps for analysis. 



A. Free energy of forming an ion-water cluster in 
the presence of the bulk medium 

We calculate ^aq("), the free energy of inserting an 
ion within a cluster in presence of the external solvent, 
for n = 1 — 6 coordination states. For the ions studied 
here, the local chemical contribution k^T Ihxq is fully 
accounted for by this coordination number^ and thus 
these are the states of most interest. 

£,aqin) is obtained indirectly by calculating the re- 
sponse of the medium, ^l^^^^{n), and subtracting it from 
jjf^in), the hydration free energy of the ion in the n*'* 
coordination state (Eq. [T5|) . Values of /^^^('t.) were taken 
from our earlier work*^. To obtain ^'^^^^(n), we simu- 
late an ion water system in which the number of solvent 
molecules within the coordination shell was held at n. 
Further, solvent molecules within the coordination sphere 
are fully-coupled with the ion, but the ion interacts with 
the bulk medium at various fractional charge states gC, 
where q is the charge on the ion and C, — 0, 1/2 ± -\/T2, 1. 
Ion-solvent pair correlation function was obtained from 
the C = 1 simulation. 

The electrostatic contribution to ^'^^^^{n) is deter- 
mined by integrating the average potential, {(t>)q(, at the 
centre of the ion due to the bulk medium. is the 
charge seen by the bulk water molecules. (We emphasize 
that the water molecules inside the coordination sphere 
always see a charge of q.) The electrostatic contribution 
to /ioutcr(^) is approximated by the two point Gauss- 
Legendre quadrature 

A^o?itcr,clcc("-) = 2^i'l')l/2+l/VT2 ~^ (<^) l/2-l/\/T2 ) 

which is exact to fourth order in perturbation theorj*^. 

For obtaining the van der Waals contribution to 
A*outcr('^)j ^6 determine the distribution of ion-bulk in- 
teraction energies from the C simulation. The van 
der Waals contribution to fI^'!|^^^^{n) is then obtained by 
approximating this distribution by a Gaussian^^ within 
the inverse form of the potential distribution theorem. 



B. Free energy of forming an ion- water cluster in 
the absence of the bulk medium 

Ion-water clusters within the coordination sphere were 
studied at different fractional charge states to obtain the 
electrostatic contribution to the free energy of forming 
an ion-water cluster when there is no bulk medium out- 
side the coordination sphere. (Ion-solvent pair correla- 
tion function was obtained from the C = 1 simulation.) 
To calculate the van der Waals contribution to the free 
energy, we first obtain the distribution of ion-solvent in- 
teraction energies for the neutral ion. We obtain the un- 
coupled binding energy distribution by performing test 
particle insertions of the neutral ion in neat water clus- 
ters. The vdW contribution is then obtained by his- 
togram overlap^. 
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IV. RESULTS AND DISCUSSION 

A. Pair correlation in K[H20]6(aq) and K[H20]6(g) 
clusters 

For hydrated ions, an implicit assumption in using 
gas-phase clusters to calculate solution phase thermody- 
namic properties is that the average distribution of water 
molecules in the solution phase and gas phase clusters are 
similar. Examining the ion- water pair correlation (Fig.[T]) 
does validate this assumption, but only for the low coor- 
dination states. 




3.0 3.5 

r (A) 



FIG. 1. The distribution, gKo{r), of water around in 
the presence (aq) and absence (g) of the external medium for 
the K[H2 0]g cluster. The radius of the coordination sphere 
A = 3.7 A. The pair correlation in presence of the medium 
shows enhanced structure relative to that for the cluster in the 
absence of the bulk. This suggests that the external medium 
promotes a better packing of solvent around K^. Inset: The 
distribution, of the farthest water molecule from K"*". 

Observe that in the absence of the bulk medium, this water 
molecule is closer to the boundary of the cluster. 



Figured] shows that the presence of the medium causes 
the cluster to be better packed, resulting in a somewhat 
sharper pair correlation for the cluster in the bulk than 
in its absence. As anticipated, the effect of the medium is 
most pronounced for the water molecule that is farthest 
from the ion (Fig. [I] Inset). 

The distribution of the water molecules comprising the 
n = 1, . . . , 3 coordination states is insensitive to the pres- 
ence of the external medium (data not shown), a fea- 
ture that is in accordance with our earlier— finding that 
there is a core group of water molecules with which the 
ion interacts strongly enough that the effect of the bulk 
medium on these water molecules is small. The higher 
coordination states are however influenced by the bulk 
medium as well. For the cluster with no bulk medium 
outside, the 6"^ water molecule is closer to the bound- 
ary of the cluster, a position that also allows this water 
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FIG. 2. The average solvent-solvent interaction energy, 
(eiu-ui), of ion- water clusters for different coordination states. 
For curves denoted by 'aq' (filled symbols), {sw-w) is obtained 
by extracting the cluster from simulations in the presence of 
the bulk medium. Curves denoted by 'g' (open symbols) cor- 
respond to clusters simulated in the absence of a bulk. The 
coordination radius A = 3.7 A. For higher coordination states, 
the average excess internal energy is lower in the case of a 
cluster without a bulk medium outside. 



molecule to associate with the remaining water molecules 
in the cluster rather than with the ion. Similar incom- 
plete shell effects, where the second hydration shell starts 
to form before the first shell is complete, have been in- 
ferred on the basis of experimental studies on gas phase 
clusters with water— or ammonia52*S2, as solvent ligands. 

Figure [5] shows that the bulk can stabilize configura- 
tions of the water molecules that would be unfavorable 
otherwise. The average excess internal energy (sw-w) of 
the n- water molecules in the coordination sphere for the 
higher coordination states, relative to clusters in the ab- 
sence of the bulk medium, is higher for the cluster that 
is extracted from simulations including the bulk. Fur- 
ther, (ew-w) is relatively insensitive to the presence of the 
bulk for the low coordination states. Thus we find that 
the medium promotes better packing of water molecules 
around the ion by stabilizing configurations of the cluster 
which are otherwise unfavorable. 



B. Medium effects on the free energy of forming 
ion-water clusters 



Wn, the free energy of forming an ion plus n- water 
molecule cluster comprises a local contribution, S,iiq{n), 
and the response of the bulk medium, /io?itcr(") ^ l^a'i^) 
(Eq. [T4|) . We next consider each of these contributions 
separately. 

FigureOshows that ^aq('^) is insensitive to the presence 
of the medium for small clusters (n = 1,2,3). However 
for larger clusters {n = 4, 5, 6), not having a bulk medium 
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FIG. 3. The local contribution, ^(n), to the free energy of 
forming an ion plus n-water molecule cluster. The filled sym- 
bols, i^aq(n-), are results in the presence of the bulk (aq) and 
the unfilled symbols, ^g(n), in its absence (g). The radius 
of the coordination sphere is 3.7 A. For n = 1,2,3, Caq(w) is 
insensitive to the presence of the bulk medium outside the co- 
ordination shell, whereas for n = 4,5, 6, (,{n) is more favorable 
with the medium than without. 



to stabilize the cluster leads to a higher free energy. The 
difference, S.a.qi'n) — ^g(n) is always negative. Thus ignor- 
ing this difference while implementing cluster variation 
will predict a less favorable contribution due to that co- 
ordination state: in effect, the probability of observing 
these higher coordination states will be predicted to be 
lower than when the difference is included. 



Our previous work showed that increasing the coordi- 
nation radius A decreases the chemical term Wn- Intu- 
itively, we expect a more favorable clustering free energy 
with increasing coordination radius because the ion can 
be better accommodated by the coordination shell sol- 
vent molecules, that is the local contribution becomes 
favorable and we expect Wn to decrease because ^aq(") 
is expected to decrease. However, we find that increasing 
the coordination radius has no effect on the local contri- 
bution, ^s^qin■), for small {n < 3) clusters and becomes 
marginally unfavorable for larger n = 4, 5, 6 clusters (Fig 
m Top). The observed decrease in Wn with increasing 
A is in fact found to arise from a increasingly favorable 
medium response: Moutcr("-) ~ Mq^(O) decreases as A in- 
creases (Fig m Bottom) . This feature emerges because 
the coordination shell solvent becomes insensitive to the 
presence of the bulk for large coordination radii. It is 
this feature that compensates for the increase in ^aq('T-) 
and leads to the observed decrease in Wn- 




A increases 



n 



FIG. 4. Top: The local contribution, 5aq('T-), of forming Na"*" 
plus n-water clusters. Bottom: The long-range contribution, 
Moutcr('^) ~ Mcf (0)i to the free energy of forming Na""" plus 
n-water clusters. Note that increasing the radius A of the 
coordination sphere does not greatly affect ^aq(n) for n < 3 
but increases it for n = 4, 5, 6. On the same scale, the response 
of the medium is more pronounced for the same change in 
A and it tends to become insensitive to the presence of the 
coordinating solvents for large A. 



C. Cluster variation and optimal coordination 
states 

Figure [S] shows the local chemical contribution (Eq. [3]) 
to the hydration free energy; a maximum term approxi- 
mation has been used. Note that Eq. [3] can be rewritten 
asi^ 

lnxo = -ln(l + ye-'5^"^), 

where Pn is the probability of observing n solvent 
molecules in the coordination sphere in the absence of 
the solute. Thus the maximum term approximation is 
simply Inxo ~ Wn — Inpn/po- When the bulk medium is 
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FIG. 5. A maximum term approximation to the local chem- 
ical contribution (Eq. lna;o « PPWn — lnp„/po- For K"*", 
neglecting the bulk medium in the local ion water clustering 
reaction leads to the identification of the n = 2, coordination 
state as the optimal (or most probable) coordination state. 
Account for the medium indicates that n = 7 is the most 
probable coordination state (cf Fig. 7 in Ref . [i^) . For Na^, 
in the presence (absence) of the medium n = 7 (n — 4) states 
are indicated to be optimal (cf. Fig. 3 in Ref |43) . 



present, £,aq{n,) defines the local contribution to Wn and 
when it is absent, ^g(n) defines the local contribution. 
In both cases, we use the same long-range contribution 

For K"'' , Fig. [5] shows that neglecting the role of the 
bulk medium in the local ion-water clustering shifts the 
predicted most probable coordination state to ti = 3, 
whereas the value obtained in simulations is n = 7 for 
A = 3.7 A: within the maximum term approximation, 
is a minimum for rt = 3 (n = 7) in the absence (pres- 
ence) of the medium. For Na"*", neglecting the medium 
suggests an optimal coordination state of n = 4, whereas 
in the presence of the bulk, a value of n = 6 is predicted. 



The value obtained in simulations^ is n = 6. For both 
Na"*" and K"*" , any discrepancy between the optimal coor- 
dination predicted using cluster variation and those ob- 
served in simulations is well within the uncertainty in the 
calculations. 

Earlier— we had defined dominant hydration struc- 
tures on the basis of how an increment in the coordination 
number contributes to fceTlnxo, the local chemical con- 
tribution to hydration. If in going from n — > n -t- 1, the 
contribution to k^T Inxo (Eq. [3]) were only of the order 
of thermal energies (and substantially smaller than the 
contribution obtained in going from n — 1 — ^ n), then 
the n-coordinate state was regarded as dominant. Since 
this definition is closely tied to the local interactions, 
it is expected to be insensitive to the definition of any 
chemically reasonable coordination radii. Fig. [S] shows 
that past n = 4 for both Na+ and K+, the change in 



kBTlnxo is only modest. On this basis we would con- 
clude that for the potential model used here, n = 4 is 
the dominant coordination state for Na+ and K+ , as was 
found earlier—. 



V. CONCLUDING DISCUSSIONS 

Effort has focused on using thermochemical data ob- 
tained either experimentally or from quantum chemical 
calculations from gas phase clusters to estimate single ion 
hydration free energies. An assumption often implicit in 
these studies is that the configurations of the clusters in 
the gas phase are similar to their aqueous counterparts. 
Our results show that this assumption is only true for 
small cluster sizes. For the Na+, K+, and F~ ions con- 
sidered here, clusters with less than three water molecules 
satisfy this requirement. For larger clusters, the external 
medium starts to influence the local ion-water interac- 
tion. The external medium stabilizes ion-water cluster 
configurations that are better packed around the ion, such 
that the local ion-solvent interaction contribution to the 
free energy of cluster formation is more favorable in the 
presence of the medium than without. For higher coor- 
dination states of the ion, and, more generally, for the 
coordination structure and thermodynamics of soft ions, 
our results show that the medium will play a sizable role 
in the coordination structure and thermodynamics of the 
hydrated ion. For these cases, accounting for the molec- 
ular characteristics of the bulk medium on the ion- water 
cluster is important in inferring the structure and ther- 
modynamics of the hydrated ion. 



Appendix A 

In Ref. |4^. the long-range contribution to Kn/Kn was 
erroneously left out from the equations. But these con- 
tributions were all correctly considered in the numerical 
work and hence no result is affected. Since presenting 
these corrections also provides a helpful alternative per- 
spective on the equilibrium constants appearing in the 
quasichemical theory, we note those corrections here. 

Eqs. 10 in Ref. |43 should be rewritten as 



In 



IV. Z, 
1 



w V 



7,1 



e-^^"+ie-^^(^"'«d7e" (Ala) 

g-/3C/„g-/30(K",/3)^^„ , (Alb) 



where 7„ = e-'^(^™tcr(")-A'r (o)) (Eq. [TS]) and is the 
configurational partition function of a water molecules. 
Un+i is the potential energy of the ion plus n solvent 
molecules within the coordination volume and C/„ is the 
potential energy of the solvent molecules in the absence 
of the ions. 0(7?.", /?) is the field of the bulk medium on 
the n solvent ligands in the coordination volume in the 
absence of the ion. Eq. lAlal includes the factor 7„ that 
was missing in Eq. 10a in Ref. 
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Based on the above definitions, we can show that 

= IZ7»<e"''"''l">oP... (A2) 

n 



where once again the factor 7„ was left out in Ref. H^. 
AJ7 = Un+i — Un is the interaction energy of the ion with 
the n solvent molecules within the coordination sphere. 
Note that ^aq(n) = -l3ln{e-'^^'^\n)o (Eq. [H]). Eq. ED 
includes the factor 7„ that was missing in Eq. 11, Ref. 

E3. 

Finally, we note that Eq. A7 in Ref. is properly 

MoutGr(^) f-Ild not /^outor 

(0) as was indicated. 
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